Quantum Expression of Classical Chaos 

S ' Sudhir R. Jain*^^\ Benoit Gremaud^^-* , Avinash Khare*-^-* 

^ ' (^'' Theoretical Physics Division 

Bhabha Atomic Research Centre 
Mumbai 400094, India 

m ; 

CN ■ (^) Laboratoire Kastler Brossel 

Universite Pierre et Marie Curie 
Q ! T12 El, 4 place Jussieu 

U ; 75252 Paris Cedex 05, France 

(^^ Institute of Physics, Sachivalaya Marg 
Bhubaneswar 751005, India 



Abstract 

VO ' Even as we understand for long that the world is quantal and buried 

^^ , in it is classical dynamics which is chaotic, finding cigcnfunctions an- 

alytically from the the Schrodinger equation has turned out to be a 
near-impossibility. Here, we discover a class of chaotic quantum sys- 
tems for which we obtain some analytically exact cigcnfunctions in 
closed form. This paves way to an exact classical and quantum me- 
chanical treatment of chaotic systems. Furthermore, we bring out con- 
r^ , nections, underlying the discovery, between different areas of physics 

and mathematics related to universality in fluctuations observed in a 
wide variety of complex quantum systems. 
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World around us is rather complex where nonlinear phenomena abound. 
Nonlinearities give birth to chaos and make it impossible to predict long- 
time dynamics (1,2). Chaotic behaviour is characterized by the existence 
of positive Lyapunov exponents, which determine the rate of exponential 
separation of very close trajectories in the phase space of the system. Upon 
casting the chaotic systems in a quantum mechanical framework, impressions 
of chaos are found in variety of statistical properties of energy levels and 
eigenfunctions (3-5). The fluctuation properties of energy level sequences of 
chaotic quantum systems agree very well with the results in random matrix 
theory (RMT) (6). These extensive studies have given place to quantum 
chaology or what we call 'statistical quantum chaos'. 



Eigenfunctions and eigenvalues contain all the information about a time- 
independent quantum system. Chaotic eigenfunctions have been studied in 
detail and there are two main ideas around which the general understand- 
ing has evolved. According to Berry (1)^ for chaotic systems, eigenfunctions 
corresponding to excited states are conjectured to be well-represented as a 
random superposition of plane waves. This conjecture has a lot of numerical 
support on one hand, and, is connected to statistical mechanics on the other 
(8,9). The other idea ensues from Heller's discovery (10,11) of scarring of 
eigenfunctions by periodic orbits, where the probability density is consider- 
ably enhanced on the periodic orbit in configuration or phase space. This 
discovery has helped in understanding how classical periodic orbits form 
the underlying fabric for quantum states, and it helps in appreciating the 
beautiful nodal patterns and contour plots of chaotic eigenfunctions. Nev- 
ertheless, an analytical expression isn't obtained, and. Berry's conjecture 
does not help in seeing how coefficients in the superposition arrange to give 
patterns. 

The objective we set out for ourselves is to get analytically exact expres- 
sions for chaotic eigenfunctions in closed form. As we shall see, the resulting 
aesthetic beauty is due to some 'magic' connections (12) involving RMT and 
its relation to other topics in physics and mathematics, including complex 
quantum systems, Riemann zeta function, exactly solvable many-body prob- 
lems, partial differential equations and the Riemann-Hilbert problem (13). 
These different areas are related to each other by the statistical properties of 
sequences characterizing them, like energy levels of nuclei and chaotic quan- 
tum systems, zeros of the Riemann zeta function, and eigenvalues of transfer 
matrices of disordered conductors. Moreover, joint probability distribution 
function (JPDF) of eigenvalues of random matrices is given by a statistical 
mechanics problem involving particles with Coulomb-like interactions (14)- 
The JPDF gives the probability with which eigenvalues {Ei\ are found in 
intervals [Ei,Ei + dEi]; all the correlation functions giving various physical 
quantities follow from this. These problems belong to the general class of 
exactly solvable statistical mechanics models. In turn, these models can be 
converted to problems in Hamiltonian dynamics where one can study the 
classical and quantal aspects. 

In this paper, we extend these connections to another class of random 
matrix models developed to explain intermediate statistics found in pseu- 
dointegrable systems, the Anderson model in three dimensions at the metal- 
insulator transition point and in certain problems in atomic physics (15-18). 
The sequences involved in these systems exhibit a behaviour intermediate to 
regular and chaotic. Once again, the eigenvalue distribution of such random 
matrices is related to an exactly solvable statistical mechanics model where 
the interaction is a screened Coulomb-like along with an additional three- 
body term (19). In this article, the abovementioned model is mapped into 
a problem in Hamiltonian dynamics in d dimensions which is shown to dis- 



play chaos. Thus, the exactly solved many-body problem is non-integrable 
in the sense that there are lesser number of constants of the motion than 
the degrees of freedom. Quantum mechanically, some chaotic eigenfunctions 
are then found analytically even though the underlying classical dynamics 
at those energies is chaotic. This is illustrated through Sections 2 to 4, the 
ensuing general perspective is given in Section 5. 

It is worth noting that the existence of Bose-Einstein condensation in a 
related one-dimensional many-body problem at zero temperature is proved 
recently (20). To prove the existence of a Bose-Einstein condensate in one 
dimension at non-zero temperatures, we need the excited states. Thus, the 
development presented here is of importance to the general theory of quan- 
tum phase transitions. 

For quantum cat maps (toral automorphisms), analytical form for the 
eigenfunctions is known (21,22), where the solutions were possible because 
the semiclassical studies turned out to be exact. These are very interesting 
results but the methods are so specialized that they do not throw much light 
on other related problems. In Section 5, we discuss the possible advances 
emerging from the discussions that follow in the context of general relevance 
of the results and the theme of this article. 

2 From Many-Body Problem to One-Body Prob- 
lem 

In 1999, a many-body problem in one dimension was discovered (19,23) 
where the nearest neighbours interact via a repulsive interaction which is 
inverse- square in distance between the particles and an attractive three- 
body interaction, also inverse-square in position coordinates of the particles. 
The A^-particle problem on a circle has the Hamiltonian : 

N 2 2 N 

- '^;^Ecot 



"vr 

z 


{Xi - 


Xi+l) 




..) 


cot 




-Xi+l) 



vr 

-{Xi^l-Xi) cot -(Xj-Xj+l) , (1) 



with Xj = XN+j- In the rest of this article, we will take the mass as unity 
and the circumference L = vr, for notational simplicity. The potential is 
singular whenever Xj+i = Xj+mr (n being an integer) with an inverse-square 
singularity. This leads to disconnected domains where the wavefunctions are 
zero at singular (hyper-)planes in the quantum problem. We choose the 
domain in which the particles are ordered as xi < X2 < ... < xi + n. The 
center of mass (CM) motion can be separated in this case by writing the 
amplitudes of the motion of particles around the CM, X = -^ J2i=i Xi in 
terms of normal mode coordinates (24)- Thus, we can write the positions as 



i=X + yj, 
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for even N {M = N/2); and for odd N, without the last term on the right 
hand side and summation going upto n = M. The quantum Hamiltonian 
operator transforms to 

H = Hcm{X) + -ffbilliard({9n}) 

1 a2 -I M o2 N 

where W^[yj({Qn})] is the potential term in (1). Note that n = ±1, ±2, ..., itM 
for odd N while n = ±1, ±2, ..., ±(M - 1),M for even N. Thus, i^biUiard 
represents a single particle in an (A^ — 1)- dimensional domain bounded by 
the (hyper-)planes where potential becomes singular forcing all the wave 
functions to be zero on the boundaries of the domain. This gives us a class 
of billiards as a function of /3 (defined through G = 0^ and g = f5{j3 — l)) and 
from dimensions two to {N — 1). A similar connection was realized by Rey 



and Choquard [24| when they considered the Calogero-Sutherland-Moser 
(CSM) system and showed that the A^-body problem in this case is mapped 
to an integrable billiard problem in A^ — 1 dimensions. It is interesting to 
note that the CSM and our model coincide for three particles as nearest 
neighbours are all the neighbours and the cotangent term is unity owing to 
a trigonometric identity. Unlike this, and most interestingly, our many-body 
problem leads to non-integrable billiards for particles greater than three 
and hence, billiards of dimension greater than two. We now turn to show 
that for A^ > 3, our billiard models are, in fact, chaotic. This also implies 
non-integrability of the many-body problem. 

3 Family of Classical Billiards 

For the sake of concreteness, we concentrate our discussion on A^ = 4, leading 
to a three-dimensional billiard. Equations of motion are singular each time a 
collision occurs (involving 2, 3 or 4 particles), but two-particle collisions are 
forbidden by energy conservation, giving rise to smooth integration using 
standard Runge-Kutta method, until one reaches a multiple collision (25). 
Preliminary studies show that regularization of the classical motion near 
these collisions could be done, but its exact implementation has not been 
achieved yet. Nevertheless, we believe that this problem does not alter results 
presented in this section. For each trajectory, we have also computed the 



associated monodromy matrix, M (whose symplectic structure is used as a 
relevant test), allowing us to extract evidence of chaotic behaviour in our 
system. For this, we have computed the "Liapunov exponents", A(T) = 
ln(|M(T) -eol) using three different vectors : eg, eg- and eg, which are, 
respectively, unit vector parallel to the flow at an initial time, unit vector 
perpendicular to the energy shell at an initial time and a "random" unit 
vector, namely l/\/6(l, 1, 1, 1, 1, 1). A"(r), being equal to In |X(r)|/|X(0)|, 
is bounded entailing thereby a vanishing Liapunov exponent, providing a 
reference scale for further numerical estimation of non- vanishing Liapunov 
exponent. For /? = 2, at classical energy equal to the quantum ground state 
energy eg = 4/?-^ (see Eq. (q)), and initial conditions, 

qi{0) \ ( 0.1000 \ / pi(0) \ / LOOOO 
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results on Liapunov exponents are plotted in Fig. ??. As expected, behavior 
of A"{T) is substantially different from those of A-'-(T) and y{T), emphasiz- 
ing thus the presence of hard chaos in the system, with Liapunov exponent, 
A ~ 6 (26). 

A more appealing evidence of hard chaos is obtained on plotting Poincare 
surface of sections (PSOS). For a generic 3D time-independent system, their 
dimensionality (4D) makes them quite useless for visualising. Fortunately, 
in our case, symmetry properties of the Hamiltonian, viz., invariance under 
(7i <-^ g_i exchange allows us to consider the reduced phase space made of 
trajectories for which qi = g_i and pi = p-i at any time, leading thus 
to an effective 2D system. Fig. ?? depicts, for /5 = 2, the reduced PSOS 
defined by q2 = 0, at different energies : Ea = eo = W (ground state energy, 
see Eq. (|)), Eh = ei = 25.5 (see Eq. (|lO|)), Ec = en = 36 (see Eq. (^) 



and Ed = 100. The empty areas appearing on all these plots correspond to 
trajectories having a four body collision in their past, this has been checked 
by propagating the corresponding initial conditions backward in time (27). 
Nevertheless the system appears to be fully chaotic in the reduced phase 
space, which is emphasized by the fact that all periodic orbits of the reduced 
dynamics are unstable (i.e., non-trivial eigenvalues of the monodromy matrix 
are not on the unit circle). Actually these orbits are also unstable when 
considered in the full phase space. Moreover, there are also unstable periodic 
orbits not belonging to the reduced phase space. Even more, we find unstable 
periodic orbits for which the four non-trivial eigenvalues of the monodromy 
matrix form a quadruplet (A, A*, A~^, A~^*), A being complex. For example, 
for /3 = 2, there is a periodic orbit of length 1.1083386 with initial conditions 
(classical energy is eo) : 

/ 9i(0) \ / -0.228018 \ / pi(0) \ / 2.06029 

q_^{0) = -0.261062 P-i(O) = 0.737421 | (5) 
y 92(0) J \ 0.000000 ) V P2(0) / \ 5.07910 



for which we find A = —367.64192 + /.64. 633456, this is shown in Fig. ??. 

/3 = 2 is not a special value, we have checked that all results presented 
in this section hold for other values /? > 1, with eventually apparition of a 
mixed dynamics for large (5. A complete study of classical dynamics has to be 
done, especially a systematic search for periodic orbits and their properties 
(variations with both (3 and energy, bifurcations, and so on). 

4 Family of Quantum Billiards 

We shall now use the known exact energy eigenstates (28) of the newly 
discovered A^-body Hamiltonian (1) and obtain few energy eigenvalues and 
eigenfunctions of the corresponding (A^— l)-dimensional billiard. The billiard 
eigenfunctions are obtained by eliminating the center-of-mass dependence 
while the eigenvalues are determined by subtracting off the center-of-mass 
energy. 

Let £'fe , ipk denote the energy eigenvalues and eigenfunctions of the A^- 
body problem (1) with periodic boundary conditions, i.e. iJVfc = E^ipk- 
Then the exact ground state is given by (19,23) 

N 

V^o = n I s™(^i - ^i+i) 1^ ' ^0 = A^/3^ (6) 

j 

provided g, G are related by p = /3(/3 — 1), G = /3^. 

As mentioned in the Introduction, -(/'q exactly gives the analytical form 
of the JPDF of eigenvalues of random matrix ensemble relevant for inter- 
mediate statistics (15,18,19) observed in plane polygonal billiards (15,29), 
Aharonov-Bohm billiards (17), Anderson model in three dimensions (16), 
and so on. 

In addition to the ground state, a few of the excited energy eigenstates 
have also been obtained in this case (28) and are given by {ipk = tpofpk) 

ct)i = ei , ^1 = So + 2 + 4/3 , 

(t)N-i = CTv-i , ^TV-i = Eq + 2N -2 + AI3 , 

N 
^N = eiCN - Y—JgeN , En = Eo + 2N + 4 + 8/3 . (7) 

Here Cj, (j = 1,2, ...,A^) is an elementary symmetric function of order j 
in the variable Zj. For example 62 = Z1Z2 + Z1Z3 + ... (containing A^(A^ — 
l)/2 terms) and Zj = exp(2iXj). Note that (pk is an eigenfunction of the 
momentum operator with eigenvalue k (30). Following the treatment in (24) 
for the CSM billiard, it is easily shown that if the billiard Hamiltonian 
satisfies the eigenvalue equation HsXh = CfcXfe) then the eigenvalues e^ and 
the eigenfunctions Xk are related to Ek,ipk by Xk = exp{—2ikX)'ipk with 



ek = Ej. — j^. Finally, we then obtain the following excited eigenvalues and 
eigenfunctions for the (A^ — l)-billiard 



Xic = ipo 
Xls = V'o 



2 
cos -r7([-^ — ^xi — X2 — ... — xn) + cyclic permutations 

2 
sin — ([-^ — ^]xi — X2 — ■■■ — Xn) + cyclic permutations 



,(8) 
,(9) 



eic = ei, = ^0 + 4/3 + 2 - - , (10) 

^ Nf3 

XN = 1pO Yl COs2{xi - Xj) + ^ , (11) 

eN = Eo + 8P + i, (12) 

where Eo is as given by Eq. (|^. Note that even though xn is A^-dependent, 
the energy difference ejsf — Eq is A^-independent. 

These expressions conclude our illustration of the connections between 
different areas of research in physics and mathematics, as noted in the In- 
troduction and by Kadanoff (12). The classical dynamics at energy equal to 
the energy eigenvalues of the eigenfunctions is chaotic, yet the eigenfunctions 
are not random superposition of plane waves. However, since /? can assume 
any value, the functions could be quite complicated, yet comprehensible. 

5 Summary 

We underline a fundamental theme which is at the heart of our findings. Let 
us first present two scenarios - Scenario 1: From pseudointegrable systems to 
a class of fully chaotic systems, and. Scenario 2: From fully chaotic systems 
to a class of integrable systems. 

Scenario 1 

Let us begin with classically pseudo-integrable billiards which are nei- 
ther chaotic (zero Kolmogorov-Sinai (KS) entropy) nor integrable. On quan- 
tization, we get the energy levels and eigenfunctions. The local statistical 
properties of energy levels agree well with the Short-Range Dyson Model 
(SRDM). The JPDF of eigenvalues is such that its analytical form coincides 
with the square of the ground-state wavefunction of the many-body problem 
(Cf. Section 2). As shown in Section 2, this many-body problem is mapped 
to a family of quantum billiards. The classical analogues of these quantum 
billiards are fully chaotic. 



Scenario 2 

We begin with classically chaotic systems and note that on quantization, 
the local fluctuation properties agree well with those of canonical RMT. 
This is very well-studied. The JPDF of eigenvalues of the random matrix 
ensembles is such that its analytical form coincides with the square of the 
ground-state wavefunction of the CSM. In the CSM, there are A^ particles 
on a circle and each pair interacts through an inverse-square potential. This 
is a completely integrable system. Since the ground-state is crystalline, here 
again the mapping to billiards in (N — 1) dimensions is possible. However, 
in this case, the billiards are fully integrable for any N. 

Thus, we see that the two random matrix models - canonical and the 
SRDM - both have two arms, one arm leads to chaos {hxs > 0) and the other 
to order {hxs = 0). The arm which leads to analytical results for chaotic 
systems belongs to the SRDM and the associated many-body problem due to 
Jain and Khare (19). Thus, we suggest that for quantizing chaotic systems, 
the route through many-body problems may be more worthwhile. Of course, 
it would require a great deal of ingenuity to construct such many-body 
problems. However, the recent work of Glashow and Mittag (31) shows that 
it is possible for any triangular billiard, where the authors generalise an old 
connection of Onsager and Sinai (32). 

There is another important future direction which suggests itself from the 
Scenarios presented above. The quantum systems that we usually deal (bil- 
liards, for instance) with, have an inflnite-dimensional Hilbert space. Thus 
the Hamiltonian matrix is infinite-dimensional. This means that the random 
matrices also should be of an infinite order, which manifests as universality 
in the scaling limit. This would lead to many-body problem with infinite 
degrees of freedom - known as fields or classical fluids governed by par- 
tial differential equations. An instance of relevance here is the Korteweg-de 
Vries equation associated with the CSM. It is only the truncations and pro- 
jections of these field solutions and the vast variety of solutions admitted by 
the partial differential equations that manifest themselves in the standard 
parlance and certainly in numerical works. As discussed by Deift (13), the 
solvability of certain partial differential equations demands the solution of 
the Riemann-Hilbert problem which, in turn, brings in matrix theory. We 
believe that herein underlie some deep secrets of physical theories. 

Finally, a concluding remark on quantum theory of chaotic systems. 
We have seen above that low-lying eigenfunctions built on classical chaos 
- termed as 'Chaotic States' - are not random superpositions of plane waves. 
Thus, Berry's conjecture would be expectedly true for highly excited chaotic 
states, and not for chaotic states in general. We believe that chaos in quan- 
tum wavefunction, even at low energies, would show up as the system evolves. 
The time-dependent wavefunction, written as a superposition of eigenfunc- 



tions, has coefficients displaying chaos. Thus, eigenfunctions form an invari- 
ant set in quantum theory, much in the same way as periodic orbits and 
fixed points do in classical theory. The most important evidence is shown 
by Heller's discovery of scars. Thus, we state - scarring of wavefunctions on 
periodic orbits is a reminder of the classical fact that on its way, an arbitrary 
trajectory is shadowed by periodic orbits. 
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Figure Captions 

1 A*(r) = ln(|M(r) -e^l), for i =\\ (continuous line), _L (dotted line) 
and r (dashed line), for (5 = 2 (see text for definition and initial condi- 
tions). As expected, A'I(T) is bounded, whereas A^(T) and \^{T) show 
linear behavior, emphasizing presence of hard chaos in the system. 



Reduced PSOS defined by g2 = 0, at different energies : Ea = eo = IQ 
(ground state energy, see Eqs. (^)), Eh = ei = 25.5 (see Eq. (p^), 
Ec = ^N = 36 (see Eq. (|l^) ) and Ed = 100. The reduced dynamics is 
clearly chaotic. (See text for discussion about areas appearing in the 
PSOS). 

Positions (top plot) and momenta (bottom plot) as function of time for 
an unstable periodic orbit of the system with /? = 2 at classical energy 
E = €0 (ground state energy). Continuous line: {qi,pi), dotted line 
(g'_i,p_i) and dashed line {q2,P2)- The four non-trivial eigenvalues of 
the monodromy matrix forms the quadruplet (A, A*, A^"*^, A^"*^*), with 
A = -367.64192 + ^64.633456 
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